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Abstract 

Recently, Duminil-Copin and Smirnov proved a long-standing conjecture by Nien- 
huis that the connective constant of self-avoiding walks on the honeycomb lattice is 
\Jl -\- \J2. A key identity used in that proof depends on the existence of a parafermionic 
observable for self-avoiding walks on the honeycomb lattice. Despite the absence of a 
corresponding observable for SAW on the square and triangular lattices, we show that 
in the limit of large lattices, some of the consequences observed on the honeycomb 
lattice persist on other lattices. This permits the accurate estimation, though not 
an exact evaluation, of certain critical amplitudes, as well as critical points, for these 
lattices. For the honeycomb lattice an exact amplitude for loops is proved. 

1 Introduction 

In 2010 Duminil-Copin and Smirnov [6] proved that the critical point of honeycomb 
lattice SAW is Zc = l/\/2 + \f2, as conjectured by Nienhuis [TJ in 1982. Their proof 
rested on establishing a connection between three generating functions for walks in a 
regular trapezoidal sublattice St,l of the honeycomb lattice, as shown in Figure [l] of 
width T and left-height 2L. All walks start on the half-edge a incident on the left wall 
labelled a. 

They first identified a parafermionic observable: 



Here a; is a point (specifically, the mid-point between two adjacent vertices) in the 
domain, ct G M, z > 0, and w is a path starting at half-edge a and finishing at point x. 
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l(io) is the length of the path and W^{a, x) is the winding, or total rotation in radians, 
when uj is traversed from a to x. 

In the special case z = Zc = 1/ \/2 + \/2 and cr ~ 5/8, F satisfies half of the discrete 
Cauchy-Riemann equations [IS] , 

{p ~ v)F{p) + {y- v)F{q) + {r- v)F{r) - 0, (1) 

where p, g, r are the mid-edges of the three edges incident on vertex v. 

By summing ([l} over all vertices in St,l, contributions from interior mid-edges 
cancel and we are left with a relation involving only mid-edges on the boundary of 
St,l- Walks from half-edge a to the boundary fall into three classes, depending on 
which boundary edge they terminate at. This observation gives rise to three generating 
functions: 



where the sums are over all self-avoiding walks from a to the a, /3 or e,e boundaries, 
respectively. 

Duminil-Copin and Smirnov showed that the relation involving boundary mid-edges 
can be written in the form 

1 ^ cos (^Y^ At,l{zc) + cos Et,l{zc) + Bt,l{zc)- (2) 
Taking the limit L — >■ oo, this becomes 

1 = cos At{zc) + cos (J) Et{zc) + Bt(2c). (3) 

This is now a relation linking three generating functions for SAW in a strip of width 
T. Duminil-Copin and Smirnov then used this to prove Nienhuis's conjecture. 
In [5] it was proved that Et{zc) = 0, so can in fact be simplified to 

1 = cos At{zc) + Bt(zc). (4) 

This modification slightly simplifies the proof of Duminil-Copin and Smirnov. 

This simpler identity, involving just two generating functions, raises the question as 
to its applicability to other two dimensional lattices, such as the square and triangular 
lattices. However, it should be remarked that the critical points are not known exactly 
for these lattices, though high precision numerical estimates are available. The proof 
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of Duminil-Copin and Smirnov identifying the critical point cannot be repeated for 
these lattices, as there is no known appropriate parafermionic observable satisfying an 
identity like ([1]). As shown by Ikhlef and Cardy [10], the dilute 0(n) model on the 
square lattice does have a parafermionic observable, and this can be used to identify 
its critical point. In the n — limit, however, that model is not the usual SAW model. 
For the triangular lattice we know of no appropriate parafermionic observable. 



T 




Figure 1: Finite patch S^^i of the honeycomb lattice. Paths run from mid-edge to mid-edge 
acquiring a weight z for each step. The paths start on the central mid-edge of the left boundary 
(shown as a). 

However, leaving aside considerations of discrete holomorphicity and the existence 
of parafermions, and motivated by a more detailed study of the width (T) dependence 
of the two generating functions At{z) and Bt{z), we investigated which features of 
the identity ^ carry over to other two-dimensional lattices. 

We did this by calculating data for the generating functions At{z) and Bt{z) in 
strips, for T < 10 on the honeycomb lattice, for T < 15 on the square lattice and for 
T < 11 on the triangular lattice. In Figure[2]we show a plot oicofi{?>TT/i)AT{z) + BT{z) 
for T e [1 . . . 10], on the honeycomb lattice, showing the intersection of all curves at 
(ze, 1), as expected from Q. 

When we repeated this calculation using the square and triangular lattice data, we 
found that there was no unique point of intersection, though the discrepancy was too 
small to be visible on a plot such as Figure [2j Instead we observed a definite monotonic 
width dependence, described in greater detail below. So we first investigated a more 
general form of (g]), allowing the "constants" to be width dependent. That is to say, 
we assumed 

1 = c^{T)At{z,) + cp{T)Bt{z,). 
We fitted successive pairs of points {At{zc), Bt{zc)) and {At+i{zc), Bt+i{zc)) in order 
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Figure 2: Plot of cos(37r/8)AT(^:) + Bt{z) for T G [I-.IO], for the honeycomb lattice, showing 
interception of all curves at (zc, 1). The insert shows a close-up or the region of interception. 



to estimate Ca(T) and cp{T), and found a weak T dependence in both "constants". 
More significantly however, we conjecture that 



lim Ca{T)/cp{T) 

T— >oo 



cos(37r/8). 



just as in the honeycomb lattice case, based on agreement to more than 5 significant 
digits for both the square and triangular lattices. In hindsight, this is perhaps not too 
surprising, as the constants multiplying the two generating functions arise from the 
winding angle of contributing graphs, and these are independent of lattice for the two 
generating functions considered, being ±7r radians for A{z) and for B{z). 

Assuming this ratio is indeed cos(37r/8) for all lattices, we re-analysed the data 
with this assumption implicit. That is to say, we fitted the data to 

c(T) = cos{3n/8) At{z,) + Bt{z,). 

We were able to estimate both the limit liniT^cx) c(T), which is of course lattice depen- 
dent, and also the nature of the T dependence. Putting all this together we estimated 
for the square lattice, 

1.024966(1 - 0.14/T^) w cos(37r/8)AT(2c) + Bt{zc). (5) 

For the triangular lattice the corresponding result is 

1.901979(1 - O.l/T^) « cos(37r/8)AT(zc) + Bt{z^). (6) 
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The leading constant is expected to be sufficiently accurate so as to restrict errors to 1 
or 2 places in the last quoted digit. The correction term is given as 0(1/T^) but that 
exponent is a guess based an a numerical estimate in the range (1.9,2.1). Finally the 
magnitude of that term is claimed to be accurate only to 10% to 20%. 

This limiting behaviour as T — cx) suggests a new numerical method for estimating 
the critical point. For the honeycomb lattice the intersection point oi cos{37t/ 8) At {z) + 
Bt{z) for any two distinct values of T uniquely determines Zc- For the square and 
triangular lattices, we instead looked at the intersection point of cos(37r/8)AT(^) + 
Bt{z) and cos{iTT /%) At +i{z) + Bt+i{z). Call this intersection point Zc{T). Then one 
expects limT^oo Zc{T) = z^- 

In this way we estimated 

= 0.3790522775 ± 0.0000000005 

for the square lattice and 

Zc = 0.240917572 ± 0.000000005 

for the triangular lattice. These estimates can be compared to the best series estimates, 
based on analysis of very long polygon series, which are Zc(sq) = 0.37905227776 [TTl 
[T^ . with uncertainty in the last digit, and Zc(tr) = 0.2409175745 [13], with similar 
uncertainty. 

Another result proved by Beaton, de Gier and Guttmann [2] is that XvaiT^oa Bt{z) = 
B{z) = for z < Zc- This means that, in this limit, the identity ^ further simplifies 
to 

1 = cos(37r/8)^(zc), (7) 

where A{z) = liniT^oo At{z). Now A(z) is the generating function for loops (half-plane 
walks which begin and end on the boundary), which is expected to behave [9] as 

A{z)^ao + ai{l-z/zcf^^^, 

so ([71) implies an exact value for the critical amplitude: 

flo = 1/ cos(37r/8). 

Correspondingly, highly accurate predictions for this amplitude for the square and 
triangular lattices follow from ([5|) and ([6|), notably 

ao(sq) w 1.024966/ cos(37r/8) = 2.678365 

and 

ao(tr) « 1.901979/ cos(37r/8) 4.970111. 

We remark that, as discussed below, we have normalised these generating functions 
differently for these two lattices, so please note details of the normalization before 
using these estimates. 

Finally, we studied the behaviour of the two generating functions as T — >■ oo. We 
found Bt{z) ~ const. /T", where a ~ 0.25. This is precisely as predicted in based 
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on conjectures of Lawler, Schramm and Werner [TS] as to the number of SAWs on 
the boundary of a domain. In 6: it was pointed out that the conjecture imphes that 
Bt{zc) should decay as T~^/^ as T goes to infinity, just as we observed. 

Similarly, we can investigate how At{Zc) = At{zc) — oq decays as T tends to 
infinity. Prom Q it follows that At{zc) also decays like T-i/* and this was observed 
numerically. 



2 Honeycomb lattice 

The original identity of Duminil-Copin and Smirnov related three distinct generating 
functions for SAW in a finite domain. Letting the length L of the domain become 
unbounded changes the domain into a strip of finite width T. As proved in [2], the 
generating function Et{z) vanishes identically in that limit, so one has an identity 
relating two generating functions: 

1 = cos(3^/8)v4t(zc) + Bt{zc). (8) 
The generating functions A^iz) and Bt(z) for T finite are rational. For example, 

Mz) = Bq{z) 



, , 2z^ (l- z^ + z^ + 3z^ ~ 4z8 + zi2) 
(1 - z^ f {I - 2z^ + z'i ^ z^) ■ 

Bi{z) = 5^ 5 '-. 



(l-z4)"(l-2z2+z4 



A^{z)^^Pr\, = where 



-75z23 + 74z2-^-61z2^ + 41z29-20z3i + z^^ + 6z3^-4z3^+z39) 

Pf{z) = 2z6(l-z2) (4_4^2_g^4^8^6_4^8^^g^io_^2zi2 + i8zi4 

-10zi6 + 3zi«-3z2"-4z22 + i0z24_i0z26 + 5228_^30) 

Q2{z) = {l-z^-z' + z^~z^f X 

(l-3z2 + 3z4-5z^ + 8z«-9zio + 7zi2_8zi4^.8^i6_5^i8^3^2o_^22) 



These generating functions have simple poles, the dominant pole being at z = Zc{T) > 
Zc{T + 1) > Zc- Furthermore, limT-s-oo Zc{T) = Zc [18] . 
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Duminil-Copin and Smirnov proved that the unique solution of 

cos{Stt/ 8) At (z) + Bt{z) ^ 1, (9) 

for any T > 0, occurs at z = Zc = 1/ \/2 + ^/2. It follows then that we could work 
backwards: given only the simple rational generating functions Ao(^) a-nd Bq[z), we 
could identify the exact value of z^ simply by seeking the solution of 

cos(37r/8)ylo(^) + Bq{z) = 1. 

If we did not already know Zc, this would be a particularly simple way to find it. 

In a further demonstration of this invariant, we show in Figure [5] a plot of 
cos(37r/8)AT(z) + Bt{z) for T g [1...10], where it can be seen that the curves 
intersect at (zc, 1), in accordance with the identity ([5]). 

Let us assume that we didn't even know the Duminil-Copin and Smirnov identity 
but rather just conjectured that some linear combination of At{zc) and Bx{zc) was 
invariant. We write this invariant as XAt{zc) + Bxizc). Then by seeking the solutions, 
for z > 0, A > of the equations 

AAo(z) + Bo{z) - AAi(z) - Bi(z) = 0, 

and 

AA2(z) + B2{z) - AAi(z) - Bi(z) = 0, 

we could discover both the invariant and the exact value of the critical point from the 
exact solutions for strips of width 0, 1 and 2 given above. 

As we show below, this suggest a way to approximate Zc for other lattices, by similar 
means. We first show that, for other lattices, an appropriate linear combination of 
At{z) and Bt{z) approaches a limit as T — >■ oo, and use this observation to estimate 
the critical point for the square and triangular lattices. 

Letting the width go to infinity, another of the generating functions, -B(z), vanishes 
[2], and the identity reduces to A{zc) = 1/ cos(37r/8). Recall that A{z) is the generating 
function for loops, whose asymptotic behaviour is believed to be [5] 

A(z) ^ao + ai(l-z/z,)3/i6. 

Thus the Duminil-Copin/Smirnov identity in the limit L — >■ c» and T — ^ c» gives us 
the exact value for the amplitude term oq = A(zc) for the honeycomb lattice. 

Next we consider the behaviour of the generating functions At{z) and Bt{z) in the 
limit T — >■ c». Denote limT-s-oo Bt{z) by B{z), with a similar definition of A{z). Recall 
that, as proved in B(z) = for z < Zc. We then wish to understand exactly how 
Bt{zc) Q as T ^ oo. 

If Bt{zc) ~ const. /T", a log-log plot of Bt{zc) against T should be linear with 
slope —a as T oo. We only have data for width T < 10, so the gradient is still 
changing slightly with T in that plot. To accommodate this, we extrapolate estimates 
of the local gradient. We define the local gradient as 
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Figure 3: Plot of the local gradient of Bt{zc) and T against l/T^ ®^. The straight line is a linear 
fit to the data in the plot. 

and plot gradB(T) against l/T"-^^, where the exponent 0.85 was chosen empirically 
to make the plot linear. The plot is shown in Figure [31 and it is manifestly clear that 
the locus extrapolates to a value of a ~ 1/4. More systematic numerical extrapolation 
techniques [8] (not detailed here) lend support to this estimate. This is precisely as 
predicted [6], based on conjectures of Lawler, Schramm and Werner [15] as to the 
number of SAWs on the boundary of a domain. This is discussed in [6], where it is 
pointed out that the conjecture implies that Bt{zc) should decay as as T goes 

to infinity, just as we observe. 

From dl]) it follows that if Bt{zc) ^ cjT^I^, then At{z^ = At{z^ — ao also decays 
as T-i/'^, and this was observed numerically by a similar plot to that described in the 
preceding paragraph. 

3 Square lattice 

As discussed above, there is no parafermionic operator that applies to the SAW model 
on the square lattice or the triangular lattice, so we can't identify the critical point for 
SAW on these lattices as did Duminil-Copin and Smirnov for honeycomb SAW. We 
might, however, expect that in the limit L oo (so that we are again considering SAW 
in a strip) there should be a similar relationship between the two generating functions 
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T 


At{Zc) 


Bt{zc) 


1 


0.684928096008073 


0.760082094484555 


2 


0.825972541624066 


0.707257323612670 


3 


0.927565166390104 


0.668934606497192 


4 


1.006072923950508 


0.639202723889591 


5 


1.069537792384553 


0.615108345881821 


6 


1.122482001562161 


0.594974760428940 


7 


1.167689112421950 


0.577763265643123 


8 


1.206987841982332 


0.562788338725227 


9 


1.241640411741764 


0.549575210877016 


10 


1.272552495675558 


0.537782341996967 


11 


1.300394615482380 


0.527156358502502 


12 


1.325676196007041 


0.517504450137522 


13 


1.348792763213512 


0.508676719252903 


14 


1.370057142972426 


0.500554481834765 


15 


1.389720731591218 


0.493042273647721 



Table 1: The values of At{z) and Bt{z) at the critical point z = Zc- 



At{z) and Bt{z), with some T-dependence that vanishes as T — > oo. 
That is to say, while the relationship 

1 = cos{3t:/&) At (z,) + Bt{zc), 

which is an identity for honeycomb lattice SAW for finite width T, cannot be expected 
to hold for the square and triangular lattices, we might expect something like 

I = Cc.{T)At{zc) + C0{T)Bt{zc) 

to hold, where the "constants" CQ,(r) and cp{T) are weakly T-dependent. 

We have computed data for the square lattice generating functions in strips of 
width T, that is, At{z) and Bt{z), for T < 15, and used our best estimate 1/zc = 
2.63815853031 [IIJ[I1] to tabulate At(zc) and Bt(zc), shown in Tabled] In ,6^ these 
generating functions for the honeycomb lattice were defined to include an extra half- 
step at the beginning of the walk and at the end of the walk. This introduces an extra 
factor of z (or, as appropriate Zc,) and we have used this definition of the generating 
functions At{z) and Bt{z) for the square lattice data. 

We then fitted successive pairs of values {At{zc)^ Bt{zc)) and (At+i(zc), Bt+i{zc)) 
for r = 1 ... 14 to 

1 = c^{T)At{zc) + C0{T)Bt{z,), 

and solved the associated linear equations for Ca{T) and cp{T), using our best estimate 
of Zc- In Figures m and [5] we show plots of values of Ca{T) against 1/T^-^^ and C/3(r) 
against 1/T°-^'^. 
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We have no basis for assuming that this is the correct form we should choose to 
extrapolate these plots; rather, the T-dependence was chosen experimentally to give 
a linear plot. Extrapolated to T = oo, we find Ca ~ 0.3734 and « 0.9756. To 
obtain more precise estimates, we extrapolated these sequences using the Bulirsch- 
Stoer algorithm [5|. This algorithm requires a parameter w which can be thought of 
as a correction-to-scaling exponent. For the purpose of the current exercise, we have 
set this parameter to 1, corresponding to an analytic correction, which is appropriate 
for the two-dimensional SAW problem [4]. Our implementation of the algorithm is 
precisely as described by Monroe [TBI, and we retained 40 digit precision throughout. 
We also applied a range of standard extrapolation algorithms to the sequences {ca(T)} 
and {c/3(T)}. These were Levin's u-transform, Brezinskii's 6 algorithm, Neville tables, 
Wynn's e algorithm and the Barber- Hamer algorithm. Descriptions of these algorithms, 
and codes for their implementation, can be found in [8]. These gave results totally 
consistent with, but less precise than, those from the Bulirsch-Stoer algorithm. 

In this way we estimated = 0.373362 ± 0.000001 and c^g = 0.975644 ± 0.000002. 
Thus the ratio Ca/cp = 0.382683(2). For the honeycomb lattice the corresponding ratio 
is cos(37r/8) = 0.3826834 . . . , which is close to, and probably equal to, the square lattice 
value. We shall see in the next section that this apparent agreement also holds for the 
triangular lattice. 



0.3736 
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Figure 4: Plot of Ca{T) against 1/T^-^^ for square lattice A walks. The straight line is a linear fit 
to the last seven data-points in the plot. 

Assuming that Ca/cp — cos(37r/8) for the square lattice, we calculated elements 
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Figure 5: Plot of Ci3{T) against I/tO-SS fpj. square lattice B walks. The straight line is a linear fit 
to the last seven data-points in the plot. 



of the sequence cos{3tt / 8) At (zc) + Bt{zc) and extrapolated these using the same ex- 
trapolation method as described above. We found the limit of the sequence to be 
1.024966 ± 0.000001, compared to a value of exactly 1 for the honeycomb lattice. Us- 
ing this estimate we plotted log(cos(37r/8)AT(zc) -I- Bt{zc) — 1.024966) against logT. 
The plot displayed slight curvature, so we plotted the local gradient. 



log 



cos(37r/8)AT(^c) + Bt{zc) - 1.024966 \ 



cos(37r/8)AT-i(2c) + Bt-i{zc) - 1.024966 



T 



T-1 



against 1/T. This extrapolated to a value in the range (1.9, 2.1), so we took the central 
value and concluded that 1.024966 — |t « cos(37r/8)^T(^c) + ^t(^c) is the asymptotic 
behaviour. Finally, extrapolating estimates of the constant ci, we estimate ci sa 0.14± 
0.02. So our final result is 

1.024966 - ^ « cos(37r/8)ylT(^c) + 5t(zc), 

which is an accurate mnemonic for square lattice strips. 

For the honeycomb lattice, it has been proved [2| that limT-^oo Bt{zc) = B{zc) = 0. 
The proof applies mutatis mutandis to the square and triangular lattices. Thus in 
the limit of infinite strip width we find 1 « 0. 3733621 A(zc), giving a prediction for 
the critical amplitude A{zc) ~ 2.678365. Current series estimates (unpublished) are 
2.66 ± 0.03, some 4 orders of magnitude less accurate than this new estimate. 
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For the honeycomb lattice the intersection point of cos{3tt/ 8) At (z) + Bt{z) for 
any two distinct values of T uniquely determines Zc- For the square and triangular 
lattices, we instead looked at the intersection point of cos{3tt/ 8) At (z) + Bt{z) and 
cos{3tt/8)At+i{z) + Bt+i{z). Call this intersection point Zc{T). Then one expects 
liniT^oo Zc{T) = Zc- We extrapolated the sequence {zc{T)} using the same Bulirsch- 
Stoer method described above, and in this way we estimated 

Zc = 0.3790522775 ± 0.0000000005. 

This estimate can be compared to the best series estimates, based on analysis of very 
long polygon series, Zc{sq) = 0.37905227776 [HJIIl], with uncertainty in the last digit. 
Thus this method is seen to be a powerful new method for estimating critical points, 
giving very good accuracy, though it doesn't rival the most powerful methods based on 
series analysis of polygon series. However it does give comparable accuracy to methods 
based on series analysis of SAW (rather than SAP). 

3.1 T-dependence of the generating functions At{zc) and Bt{zc)- 

As for the honeycomb lattice, we expect Bt{zc) ^ const. /T^/"*. In Figure [H] we have 
plotted estimates of the exponent 

against l/T^ *^. This local gradient should approach —1/4 and from the plot is seen 
to do so. 

As for the honeycomb lattice, from ([8]) it follows that if Bt{zc) ^ const. /T^/^, then 
At{zc) = AT{zc) — ao also decays as T~^^'^. This was observed numerically by a similar 
plot to that described in the preceding paragraph. 

3.2 Alternative estimate of the critical point 

We showed in Section [2] that the critical point could be identified just from knowledge 
of the invariant ([5]) and the formulae for ^o(^) and Bo{z). In the case of the square 
lattice, we expect the simultaneous solution of the pair of equations 

XAr-iiz) + Bt-^i{z) - \At{z) - Bt{z) = 0, 

and 

XAt{z) + Bt{z) - \At+i{z) - Bt+i{z) - 0, 

to give a sequence of estimates of Zc{T) that should converge to the critical point z^- 
Similarly, the parameter A should converge to the ratio Ca/cfi = cos(37r/8). The merit 
of this method of estimating the critical point is that it makes no assumption about 
the value of A, while the method of estimating the critical point described in Section |3] 
assumes that A = cos(37r/8). 

We solved these equations by seeking the solution of 

{At-i{z) - At{z)){Bt+i{z) - Bt{z)) = {At{z) - At+i{z)){Bt{z) - Bt-i{z)), 
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Figure 6: Plot of gradB(T) against l/T"-^^ for square lattice B walks. The straight line is a linear 
fit to the last seven data-points in the plot. Linear extrapolation to -0.25 is well supported. 



which we call Zc{T), and then found A(r) by back substitution. The results are shown 
in Tabled 

We plotted (not shown) the estimates of Zc{T), against various powers of 1/T, 
and found a linear plot if we plotted against We extrapolated the estimates 

Zc{T) for steadily increasing T values using the same Bulirsch-Stoer extrapolation 
method described above. Rapid convergence was observed, and we estimate Zc = 
0.37905228 ± 0.00000001. This is consistent with the limit found from our previous 
method described above, though not quite as precise. 

We have similarly extrapolated the estimates of A(r), and find A « 0.38268, com- 
pared to the expected value cos(37r/8) — 0.382682. 



4 Triangular lattice 

We have also generated data for the triangular lattice in strips of widths up to and 
including 11. Using the best estimate [M] of the critical point, Zc = 0.2409175745, 
we show, in Table |31 the values of At{zc) and Bt{zc) for each strip width. For the 
triangular lattice there are two edges incident upon the origin in a strip geometry 
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T 




A(T) 


2 


0.3792132510996564 


0.3680865016160631 


3 


0.3791354275802486 


0.3724957043646600 


4 


0.3791014587212902 


0.3750327779790171 


5 


0.3790837649841775 


0.3766921607963391 


6 


0.3790736177161640 


0.3778473665510655 


7 


0.3790673896037665 


0.3786868392606266 


8 


0.3790633602596354 


0.3793176229093515 


9 


0.3790606406190476 


0.3798046222881576 


10 


0.3790587398862656 


0.3801891415440993 


11 


0.3790573721782793 


0.3804985249279112 


12 


0.3790563633515162 


0.3807514859621669 


13 


0.3790556032334455 


0.3809611990869139 


14 


0.379055019822686 


0.3811371681858631 



Table 2: Estimates of Zc{T) and A(T) for the square lattice. 



and this comphcates matters. To simplify things, we start and finish our SAW on 
the boundary in the case of the triangular lattice, in order to avoid the complications 
that arise when including an incident edge. So the extra factor of Zc included in the 
definition of these amplitudes for the square and honeycomb lattice data is not present 
in the triangular lattice data. 

As in the analysis of the square lattice data, we fitted successive pairs of values 
(At(zc), Bt(zc)) and (At+i(zc), Bt+i{zc)) for T 1 ... 10 to 

1 = Ca{T)AT{z,) + cp{T)Bt{zc), 

and solved the associated linear equations for Ca{T) and cp{T). 

To obtain precise estimates, we again applied the Bulirsch-Stoer extrapolation al- 
gorithm to the sequences {cq(T)} and {c^(T)}. Combining the results from these 
different algorithms, we estimate c„ = 0.2012028(3) and = 0.525770(3). Thus 
the ratio Ca/cp — 0.382682(3). For the honeycomb lattice the corresponding ratio 
is cos(37r/8) — 0.3826834 . . . , which (as we also saw for the square lattice) is close to, 
and probably equal to, the triangular lattice value. 

Assuming the ratio c^/cp — cos(37r/8) for the triangular lattice too, we extrapo- 
lated cos{3tt/ 8) At {zc) + Bt{zc) for increasing values of T, using our standard suite of 
extrapolation algorithms and the Bulirsch-Stoer algorithm. We estimated the limit to 
be 1.901979 ± 0.000001. We then repeated the analysis described above for the square 
lattice data mutatis mutandis, and found 

1.901979 - ^ ~ cos(37r/8)AT(zc) + St(zc). 
As remarked above, it has been proved [2] that limT->.oo Bt{zc) — B{zc) — 0. Thus in 
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T 




Bt{zc) 


1 


1 


139480549210468 


1 


457161363236105 


2 


1 


435344242350752 


1 


348134252648887 


3 


1 


641756149326264 


1 


270897362392145 


4 


1 


798515045521241 


1 


211810836367619 


5 


1 


923848231267622 


1 


164374555192450 


6 


2 


027608945103857 


1 


125001488941636 


7 


2 


115709764900265 


1 


091512525007183 


8 


2 


191966367371986 


1 


062490013670246 


9 


2 


258977760090717 


1 


036962918106255 


10 


2 


318589791981952 


1 


014238779515961 


11 


2 


372157936598986 





993807536013206 



Table 3: Estimates of At{zc) and Bt{zc) for various strip widths for the triangular lattice. 

the limit of infinite strip width we find 1.901979 « cos(37r/8)A(zj,), a prediction for the 
critical amplitude A{zc) ~ 4.970111. 

As for the square lattice case, we estimated the critical point Zc by extrapolating 
the intersection point of cos(37r/8)AT(^) + Bt{z) and cos{37r /8)At+i{z) + Bt+i{z), 
called Zc{T). One expects limT_^.oo Zc(T) — Zc- 

In this way we estimated 

= 0.240917572 ± 0.000000005 

for the triangular lattice. This can be compared to the best series estimate, based on 
analysis of very long polygon series Zc(tr) — 0.2409175745 [M] with uncertainty in the 
last quoted digit. 

4.1 T-dependence of the generating functions At{zc) and Bt{zc). 

As for the honeycomb and square lattices, we expect Bt{zc) ^ const. /T^^^. We plotted 
estimates of the exponent 

against l/T^-^^, which should approach —1/4, and were seen to do so. The figure was 
visually indistinguishable from the corresponding Figure [6] for the square lattice, so is 
not shown. 

Similarly, it follows from dH) that At{zc) = At(zc) - ao ^ At{zc) - 4.97011 decays 
as l/T^/^ as T tends to infinity. As we did for the square lattice case, we also confirmed 
this numerically. 
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4.2 Alternative estimate of the critical point 

In the previous section, we showed that, for the square lattice data, the simuhaneous 
solution of the pair of equations 

XAt-i{z) + Bt-i{z) - XAt{z) - Bt{z) = 0, 

and 

\At{z) + Bt{z) - \At+i{z) - Bt+i{z) = 0, 

gives a sequence of estimates of Zc{T) that converges to the critical point Zc- Similarly, 
the parameter A converges to the ratio Ca/cp = cos(37r/8), where the equality is con- 
jectural. We solved these equations using the triangular lattice data, by seeking the 
solution of 

{At-i{z) - At{z)){Bt+i{z) - Bt{z)) = {At{z) - At+i{z)){Bt{z) - Bt-i{z)), 

called Zc{T), and then found X{T) by back substitution. The results are shown in 
Table HI 

We plotted (not shown) the estimates of Zc{T), against various powers of 1/T, 
and found a linear plot if we plotted against We analysed the sequences in 

precisely the same way as for the corresponding square lattice data, using the Bulirsch- 
Stoer algorithm. The merit of this method of estimating the critical point is that 
it makes no assumption about the value of A, while the method of estimating the 
critical point described in Section |4] assumes that A = cos(37r/8). For the critical point 
we estimate Zc = 0.240917575 ± 0.000000005. This is of comparable precision to our 
estimate given in Section 21 but slightly less precise than the best series estimate [14] 
of Zc = 0.2409175745, with uncertainty in the last digit. 

Thus this method is again seen to be a powerful one for estimating critical points, 
giving very good accuracy. We have similarly extrapolated the estimates of A(T), and 
find A « 0.38268, compared to the expected value cos(37r/8) = 0.382682, exactly as for 
the square lattice. 



5 Enumeration of self-avoiding walks 

The algorithm we use to enumerate SAWs on the square lattice builds on the pioneering 
work of Enting [7] who enumerated square lattice self-avoiding polygons using the finite 
lattice method. More specifically our algorithm is based in large part on the one devised 
by Conway, Enting and Guttmann [5] for the enumeration of SAWs. The details of 
our algorithm can be found in [13]. Below we shall only briefly outline the basics of 
the algorithm and describe the changes made for the particular problem studied in this 
work. 

The generating function for a rectangle was calculated using transfer matrix (TM) 
techniques. The most efficient implementation of the TM algorithm generally involves 
bisecting the finite lattice with a boundary (this is just a line in the case of rectangles) 
and moving the boundary in such a way as to build up the lattice vertex by vertex as 
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T 




A(T) 


2 





241168440165255 





356318356471223 


3 





241030169141752 





366143831978748 


4 





240977832351101 





371147184391665 


5 





240953612190006 





374091365359823 


6 





240940839933527 





375992279128027 


7 





240933460889859 





377300697288292 


8 





240928899289076 





378244599187661 


9 





240925927855959 





378950553554884 


10 





240923909640445 





379493901730187 



Table 4: Estimates of Zc{T) and A(T) for the triangular lattice. 



illustrated in Figure [71 If we draw a SAW and then cut it by a line we observe that 
the partial SAW to the left of this line consists of a number of loops connecting two 
edges (we shall refer to these as loop ends) in the intersection, and pieces which are 
connected to only one edge (we call these free ends). The other end of the free piece is 
an end point of the SAW so there are at most two free ends. 

Each end of a loop is assigned one of two labels depending on whether it is the 
lower end or the upper end of a loop. Each configuration along the boundary line can 
thus be represented by a set of edge states {ci}, where 

{0 empty edge, 
1 lower loop-end, 
2 upper loop-end. 
3 free end. 

If we read from the bottom to the top, the configuration or signature S along the 
intersection of the partial SAW in Figure[7]is S = {031212120}. Since crossings aren't 
permitted this encoding uniquely describes which loop ends are connected. 

A few changes to the algorithm described in [13] are required in order to enumerate 
the restricted SAWs we study here. Most importantly the SAW must have a free end 
at the middle vertex of the top side of the strip. This is easily ensured by restricting 
the updating rules at this vertex (also signatures prior to passing this vertex can have 
at most one free end) . Specifically the middle vertex is reached when the TM boundary 
has been moved halfway through the strip. At this point the incoming edge to the left 
of the middle vertex is either empty, an upper loop-end or free. In the empty case we 
have to insert a new free end (along either the horizontal or vertical outgoing edge). In 
the upper case the loop-end is terminated and the matching lower loop-end becomes a 
free end. 

While in the free case the end is again terminated and all the edges connected to this 
free end form a SAW. However this is only a valid configuration if all other edges are 
empty since otherwise we would form configurations with more than one component. 
Secondly, in enumerating SAWs of type A the second free end must lie in the top side 
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Figure 7: A snapshot of the boundary hne (dashed hne) during the transfer matrix calculation 
of type A configurations on a strip of size 7 x 10. SAWs are enumerated by successive moves of 
the kink in the boundary line, as exemplified by the position given by the dotted line, so that one 
vertex and two edges at a time are added to the strip. To the left of the boundary line we have 
drawn an example of a partially completed SAW. The heavy lines at the top are the incoming and 
outgoing edges of the SAW. 

of the rectangle; we chose to force the free end to lie to the left of the middle vertex 
and use symmetry to count all possible configurations. In counting bridges or SAWs 
of type B the second free end must lie at the bottom of the strip. Thirdly, in [13] the 
SAWs were forced to span the rectangle, that is touch all sides, but this restriction is 
lifted in this study. 

The sum over all contributing graphs is calculated as the boundary is moved through 
the lattice. For each configuration of occupied or empty edges along the intersection 
we maintain a generating function Gs for partial walks with signature S. In exact 
enumeration studies Gs would be a truncated two- variable polynomial Gs{z) where z 
is conjugate to the number of steps. 

In a TM update each source signature S (before the boundary is moved) gives rise 
to a few new target signatures S' (after the move of the boundary line) and n = 0, 1 or 
2 new edges are inserted leading to the update Gs' {z) = Gs' {z) + z"Gs{z). Once a sig- 
nature S has been processed it can be discarded. In most studies the calculations were 
done using integer arithmetic modulo several primes with the full integer coefficients 
reconstructed at the end using the Chinese remainder theorem. Here we are not really 
interested in the exact coefficients. This makes life a little easier for us since we can use 
real coefficients with the generating functions truncated at some maximal degree M . 
The calculations were carried out using quadruple (or 128-bit) fioating-point precision 
(achieved in FORTRAN with the REAL(KIND=16) type declaration). 

In our calculations we truncated At{z) and Bt{z) at degree M = 1000 and used 
strips of half-length L — M . These choices of M and L more than suffice to ensure 
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that numerical errors are negligible as evidenced by the fact that when we solve (|3]) 
(with Za replaced by z) to find Zc for the honeycomb lattice the estimate for Zc agrees 
with the exact value to at least 30 digits, that is, to within the numerical accuracy of 
the floating-point computation itself. 

The computational complexity of the calculation required to obtain the number of 
walks of in a strip of width T and length L can be easily estimated. Time (and memory) 
requirements are basically proportional to a polynomial in M and L times the maximal 
number of signatures, iVconf, generated during the calculation. It is well established 
[12] that A^conf oc 3"^ so the algorithm has exponential computational complexity. 

The transfer-matrix algorithm is eminently suited to parallel computations and 
here we used the approach first described in [T^] and refer the interested reader to 
this publication for further detail. The bulk of the calculations for this paper were 
performed on the cluster of the NCI National Facility, which provides a peak computing 
facility to researchers in Australia. The NCI peak facility is a Sun Constellation Cluster 
with 1492 nodes in Sun X6275 blades, each containing two quad-core 2.93GHz Intel 
Nehalem CPUs with most nodes having 3GB of memory per core (24GB per node). It 
took a total of about 1800 CPU hours to calculate At{z) for T up to 15. So, the bulk 
of the time (almost 1250 hours) was spent calculating ^15(2). In this case we used 48 
processors and the split between actual calculations and communications was roughly 
2 to 1 (with quite a bit of variation from processor to processor). Smaller widths can 
be done more efficiently in that communication needs are fewer and hence not as much 
time is used for this task. 

On a technical issue we note that quad precision is not a supported data type in 
the MPI standard. So in order to pass messages containing the generating functions 
we used the MPI data type MPI-BYTE with each coefficient then having a length of 
16 bytes. 

The algorithm used for the triangular lattice is quite similar. The triangular lattice 
is represented as a square lattice with additional edges along one of the main diagonals. 
This poses an immediate problem since a boundary line drawn as in Figure [7] would 
intersect 2T edges thus greatly increasing the number of possible signatures. In this 
case it is more efficient to draw the boundary line through the vertices of the lattice. 
We then again have T intersections, however a vertex may be in an additional state 
since a partial SAW can touch the boundary line without crossing it (see [13] for further 
details). The upshot is that the computational complexity grows exponentially as 4^. 

6 Conclusion 

We started this study in order to consider to what extent, in some ill-defined sense, does 
there exist a parafermionic operator applicable to SAWs on the square and triangular 
lattice. For honeycomb lattice SAWs, ([S]) is an identity. For the square and triangular 
lattices we found, experimentally, that for walks in a strip of width T on the square 
lattice, 

1.024966(1 - 0.14/r2) w cos(37r/8)AT(zc) + Bt{z<,), (11) 
while for the triangular lattice the corresponding result is 

1.901979(1 - O.l/T^) « cos(37r/8)AT(zc) + Bt{z^). (12) 
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Since Bt{zc) — > as T — )■ oo, it follows that in the same limit At{zc) — >■ clq, 
where the numerical value of oq — l/cos(37r/8) for the honeycomb lattice, ao « 
1.024966/ cos(37r/8) for the square lattice, and oq « 1.901979/ cos(37r/8) for the tri- 
angular lattice. We provided numerical support for the conjecture that Bt{zc) ~ 
const. /T^/"* as T ^ oo, and hence that At{zc) — oq ^ const. /T^/"* also. Finally we 
show how the existence of equations ([S]), pT|) and (IT^ suggests a powerful numeri- 
cal method to estimate the critical point. In that way we found Zc exactly for the 
honeycomb lattice, and estimated Zc = 0.3790522775(5) for the square lattice and 
Zc = 0.240917573(5) for the triangular lattice. 

It has been pointed out to the authors by John Cardy (private communication) that 
in the scaling limit, all two dimensional SAW models are given by the same conformal 
field theory (CFT). Since it is known that for one of these models (i.e. honeycomb 
lattice SAW) the critical point can be found by requiring certain contour integrals to 
vanish (i.e. when summing ([T]) over a region of the lattice), it follows that in the scaling 
limit the same must be true for all two dimensional SAW. This is entirely consistent 
with our observations and the relations (fTTj) and ([T^ . 
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